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ABSTRACT 

We explore the nature of possible electromagnetic counterparts of coalescences of 
massive black hole binaries at cosmological distances detectable by the Laser Inter- 
ferometer Space Antenna (LISA). An electromagnetic precursor, during the last year 
of gravitational wave (GW)-driven inspiral, or an afterglow within few years after 
coalescence, may highlight the position in the sky of galaxies hosting LISA sources. 
We show that observations of precursors and afterglows are mutually exclusive, de- 
pending on the mass of the primary black hole. Precursors are expected to occur in 
binaries where the primary (more massive) black hole is heavier than ~ 10 7 M . They 
may correspond to on-off states of accretion, i.e., to a bright X-ray source decaying 
into quiescence before black hole coalescence, and are likely associated to disturbed 
galaxies showing signs of ongoing starbursts. Coalescences of lighter binaries, with 
masses <5x 10 6 M Q , lack of any precursor, as gas is expected to be consumed long 
before the GW-driven orbital decay. Such events would not be hosted by (massive) 
galaxies with an associated starburst, given the slow binary inspiral time compared to 
the typical time scale of starbursts. By contrast, coalescence, for such light binaries, is 
followed by an electromagnetic afterglow, i.e., an off-on accretion state rising in < 20 
yrs. Using a cosmological merger tree algorithm, we show that future X-ray missions 
such as XEUS will be able to identify, in 20 yrs operation, almost all the massive BH 
binary detectable by LISA, and, in only 5 yrs, all the LISA sources at z > 6. 

Key words: accretion, accretion discs - black hole physics - gravitational waves - 
quasars: general - galaxies: starburst 



1 INTRODUCTION 



Massive black hole (BH) binaries are considered primary 
sources of gravitational waves (GWs) detectable, through- 
out the entire Universe (up to redshifts z > 10, Hughes 2002, 
Bender & Hils 2003), by the space-based Laser Interferom- 
eter Space Antenna (LISA; Bender et al. 1994, Hils & Ben- 
der 1995), designed to operate in the low-frequency range 
(3 x 10~ 5 - 0.1 Hz). LISA will observe the signal, emitted 
during the last year of inspiral and at coalescence, by BH 
binaries of characteristic masses between lO 4 M0/(l + z) and 
10 9 M Q /(l + z) (Haehnelt 1994; Jaffe & Backer 2003; Wyithe 
& Loeb 2003; Sesana et al. 2004, 2005; Alcubierre et al. 
2005). In standard hierarchical cosmogonies, a large number 
of BH binaries is predicted to form along the cosmic his- 
tory, through collisions of (pre-) galactic structures (Begel- 
man, Blandford, & Rees 1980; Kauffmann & Haehnelt 2000; 
Volonteri, Haardt, & Madau 2003). 



LISA will operate as an all-sky monitor, and its data 
stream will record the signals from a large number of sources 
belonging to different populations of both galactic and cos- 
mological origin (Nelemans, Yungelson & Portegies Zwart 
2001; Benacquista et al. 2003; Farmer & Phinney 2003; 
Barak & Cultler 2004). Given the difficulty in disentan- 
gling the different signals, and the large number of noise 
sources, the degree of precision for the determination of 
the source position on the sky is still an open issue (Cut- 
ler 1998; Hughes 2002; Vecchio 2004; Kocsis et al. 2005). 
Analysis of LISA data alone will provide a measure of a 
number of binary parameters not corrected for the source 
redshift (Vecchio 2004; Holtz & Hughes 2005, Baker et al. 
2006). Thus, pinpointing the GW source to its electromag- 
netic (EM) counterpart within the LISA error cube would 
be of great importance. Although the detection of a GW sig- 
nal from coalescing BHs would be an extraordinary event by 
itself, the identification of an EM counterpart would greatly 
increase the scientific payoff, allowing: (1) to improve our un- 
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derstanding of the nature of the galaxy hosting coalescing 
massive BH binaries (e.g., galaxy type, colors, morphology, 
etc.), (2) to reconstruct the dynamics of the merging galax- 
ies and of their BHs, and, most importantly, (3) to break 
the degeneracy between the luminosity distance, which is 
directly measurable from the GW signal and the source red- 
shift. In this last case, LISA sources could be used (1) as 
"standard candles" ("standard sirens", Schutz 1986) to es- 
timate fundamental cosmological parameters (Schutz 2002; 
Holz & Hughes 2005), (2) to distinguish between the mea- 
sured redshifted BH masses and their rest-frame values, and 
(3), to improve our understanding of the accretion history 
of massive BHs comparing the EM luminosity with the BH 
masses and spins. 

Different possible EM counterparts of coalescing mas- 
sive BHs have been proposed in the last few years. Kocsis 
et al. (2005) proposed that LISA events are associated to 
the brightest quasar in the LISA error box, provided that 
the volume explored is limited to z < 1 (z < 3 in case of 
rapidly spinning BHs), guaranteeing an unambiguous iden- 
tification. Armitage & Natarajan (2002) suggested that a 
strong, potentially observable accretion episode should an- 
ticipate the BH collision: gas trapped inside the binary or- 
bit may produce a surface density spike, and a bright flare 
with outflows just prior coalescence. Finally, Milosavljevic 
& Phinney (2005) argued that a near-Eddington X-ray af- 
terglow should appear, delayed by few years from the LISA 
detection. On such timescale, residual circumbinary gas ac- 
cumulated during the viscous-migration of the less massive 
BH has time to refill the accretion volume, and feed the 
central BH. 

In this paper, we explore the proposed EM counter- 
parts to LISA events, relating their characteristics to the 
physical properties of binary BHs and of their host galaxies. 
We restrict our analysis to the BH binaries with mass ra- 
tios > 0.01, i.e. covering the expected range in hierarchical 
scenario of supermassive BH assembly (Sesana et al. 2005). 
Although coalescences of binaries with more extreme mass 
ratio could be important sources of GW emission, the physi- 
cal conditions for these events might be really different with 
respect to the cases here explored. We address a number of 
questions: (1) What are the BH masses for which we can 
potentially see an EM precursor or/and an EM afterglow? 
(2) Does the host galaxy show visible signs of a collisional 
interaction (a starburst) during the last-year of BH inspi- 
ral? (3) Will the next-generation of X-ray missions be able 
to identify LISA events at cosmological distances? 

The paper is organized as follows. In Sect. 2 we in- 
troduce characteristic radii and times associated to the cir- 
cumbinary disc surrounding the BHs, and describe the lat- 
est stages of orbital decay, before GW-driven coalescence. 
In Section 3 we infer relations that are necessary (not suf- 
ficient) for the occurrence of an EM precursor phase, and 
explore the possibility that tidal disruption of bound stars 
may trigger an episode of mass transfer leading to an X- 
ray flare. We then study the detectability of EM afterglows 
with XEUS. In Section 4, we discuss on the potential link be- 
tween starburst-host-galaxies and LISA events. In Section 
5 we derive our conclusions with a sketch of the domains 
where EM preglows and afterglows can occur, in the BH 
mass parameter space. 



2 BH DYNAMICS 

Coalescences of massive BHs are thought to be events as- 
sociated with mergers of (sub-) galactic structures at high 
redshifts, so that the orbital decay, from the large distance 
scale of a merger (~ 100 kpc) to the subparsec-scale, likely 
occurs in a gas-rich, dissipative environment. In cosmologi- 
cal models of structure formation the typical mass ratio of 
close BH binaries observable by LISA is ~ 0.1, with a tail 
extending down to ~ 0.01 (Sesana et al. 2005). Thus, we 
will not consider here more extreme mass ratios. 

Recently, Kazantzidis et al. (2005) explored the effect of 
gaseous dissipation in mergers between gas-rich disc galaxies 
with central BHs, using high resolution N-Body/SPH simu- 
lations. The authors show how the presence of a cool gaseous 
component is essential in order to bring, in minor mergers, 
the BHs to parsec scale separations. Gas infall deepens the 
potential well, preserving the less massive galaxy against 
tidal disruption, and leading to formation of a close BH 
pair. * Moreover, the interplay between strong gas inflows 
and star formation seems to lead naturally to the formation, 
around the two BHs, of a massive circumnuclear gaseous 
disc on scales of < 100 pc, close to the numerical resolution 
limit. Despite these advances, it is still difficult to establish 
the dynamical properties of such self-gravitating discs at 
sub-parsec scales, as well as to asses the characteristics of 
the BH orbits. t 

Escala et al. (2005; see also Escala et al. 2004) and 
Dotti, Colpi, & Haardt (2006) have recently studied the 
dynamics of double massive BHs (in the range 10 Mq < 
Mbh < 2.5 x 10 9 M Q ) orbiting inside a circumnuclear disc 
modeled with a Mestel profile. They found that orbital an- 
gular momentum losses via dynamical friction against the 
gaseous background are sufficient to pair the BHs into a 
bound Keplerian orbit (on distances of a few pc), and to 
drive initially eccentric orbits into circular before the two 
BHs form a bound state (Dotti et al. 2006). Escala et al. 
(2005) have also shown that gravitational torques excited 
by ellipsoidal deformations of disc inner regions can bring 
the binary to even closer distances. Torques can reduce the 
separation down to ~ 0.1 pc, the scale (for M > 10 8 M Q ) 
at which GW emission can bring the two BHs to the coales- 
cence in less than an Hubble time. In extrapolating the re- 
sults of the highest resolution simulations available to date, 
Escala et al. (2005) suggested that gas-driven orbital de- 
cay leads the binary to coalescence in < 10 7 yr. All these 
simulations lack of enough spatial resolution to detail the 
transition boundary separating the self-gravitating region 
of the disc from the region dominated by the gravity of the 
binary. Inside this transition zone, the disc surrounding the 
binary is expected to be accreting through viscous stress. 
If the mass ratio q (q = M2/M1 < 1) of the binary BHs 
is <C 1, then we may estimate the position of the radius 

* In collisionlcss mergers, the bulge of the lighter galaxy is tidally 
disrupted before the merger is completed, leaving its central BH 
wandering in the outskirts of the remnant galaxy (Kazantzidis ct 
al. 2005). 

t See Mayer ct al. 2006 for new preliminary simulations at higher 
resolution, which are designed to solve for the spatial and veloc- 
ity structure of the circumnuclear discs and dynamics of the BH 
pairing process. 
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7?tr at which such transition occurs, assuming that the disc 
is described by a Shakura-Sunyaev gas-pressure dominated 
Q-disc (Shakura & Sunyaev 1973), and that the disc mass 
within i? tr is of order of the binary mass. 

The disc mass in an a-disc scales with radius as 

M disc (R) ~ 3.5 x 10 8 M (^-j f E /5 R 7/5 , (1) 

where f E is the luminosity in Eddington units (computed 
using a radiative efficiency e = 0.1), and M7 is the reference 
mass of the central, more massive BH in units of 10 7 Mq. 
Thus, 

fltr ~ 8 x 10- 2 pc (1 + g) - 71 M 7 - 14 a - 57 /s°' 43 . (2) 

Given our estimate of R tI , one can fairly assume that dy- 
namical friction can shrink the BH binary down to this sep- 
aration. Indeed, on the bases of known results of disc-planet 
interactions (see e.g. Papaloizou et al. 2006 and references 
therein) we expect that M2 perturbs gravitationally the sur- 
rounding viscous gas creating a low density, hollow region, 
called "gap", which separates the outer-circumbinary disc 
from an inner gas-poor region around the primary BH. At 
a gap the gravitational torque from the secondary BH bal- 
ances the gas-viscous torque. 

In the case considered here, the appearance of the gap 
would occur at a BH binary separation 

fl gap ~ fl tr ~ 0.6i?tr, (3) 

where the prefactor 0.6 is taken from the numerical simula- 
tions of Artymowicz & Lubow (1994). The scale a tr defines 
naturally the time for binary coalescence set by the GW 
back-reaction, which corresponds to the longest time in ab- 
sence of other processes of orbital angular momentum loss: 

tGw(otr) ~ 2.5 x 10 12 yr 

xMf^V- 29 /;- 1 - 71 ?"^! + q^FieyW 

where F(e) takes into account the dependence of £gw on 
the binary eccentricity e (Peters 1964; following the re- 
sult of Dotti et al. 2006 on BH circularization, we consider 
F(e) ~ 1 hereafter). Note that at a tr , iGw < 10 Gyr for 
M > 3 x 10 7 M Q (q = 0.1, a = 0.1, f E = 1). To accomplish 
coalescence on a shorter time, additional angular momentum 
losses are necessary. 

For q "C 1, Armitage & Natarajan (2002) have shown 
that, once the gap opens, gravitational torques exerted by 
the inner edge of the circumbinary disc can still shrink the 
orbit, so that the inspiral time is set by the disc-driven mi- 
gration timescale. For their choice of binary and disc param- 
eters, coalescence is reached in ~ 10 7 yr. If we assume in our 
study that this process is efficient even for q < 0.1, and that 
migration continues down to the binary separation at which 
the disc mass inside the orbit equals M2 (Armitage private 
communication), i.e., 

a cr it ~ atr[g/(l + q)]"' 71 

, » „ 1n -2 0.71,^0.14 0.57 r-0.43 /,-n 

~ 4.7 x 10 pc q M 7 a f E , (5) 
then the time needed to migrate from a gap to a cr i t would be 

1 . w 1n 8 -0.8 ji r-0.2 f 0.4/ 1.4 1.4 \ / a \ 

torque ~ 4.2 x 10 yr a M 7 f E (a gap - a crit ) (6) 

(see eq. 1 of Armitage & Natarajan 2002). From a cr i t in- 
wards, GWs can drive the binary to coalescence on a 
timescale 



tcw(acrit) = iGw(atr)(«crit/atr) • (7) 

Note that now t G w < 10 Gyr for Mi > 2 x 10 6 Mq (q = 0.1, 
q = 0.1, f E = 1). According to eqs. 6 and 7, the coalescence 
time of the binary BHs can be approximated as 

tcoal = mhl{£Gw(atr), [^torque + *GW (cicrit )] } • (8) 



3 NUCLEAR ACTIVITY AS EM 
COUNTERPART 

3.1 Circumbinary disc accretion: the precursor 

A galaxy showing nuclear activity can be a peculiar EM 
counterpart of a LISA event. In order to have an AGN pre- 
cursor, part of the gas has to remain bound forming a small 
accretion disc at least around the more massive BH when 
the gap opens. The time of gas consumption, iduty, is given 
by 

tduty.l = Y~~^ ^dd/fi 1 M 1 + (A4disc,l(a g ap)/Ml)], (9) 

where e is the radiative efficiency, and Mdisc.i is the mass of 
the disc when the gap opens. Mdi sc ,i is obtained from eq. 1, 
using the estimate of the disc radius around Mi as given by 
Artymowicz & Lubow (1994), viz. ~ 0.45a gap for our choice 
of parameters. The AGN is still in an "on" state at binary 
coalescence only if 

Wty(a gap ) > ^coal- (10) 

The presence of an AGN precursor varies with Mi and 
q as shown in Figure 1 (for q = 0.1) and Figure 2 (for 
q — 0.01). The solid curve, labeled with f coa i = t duty in 
both figures, defines the critical BH binary separation be- 
low which gap opening still ensures nuclear activity at co- 
alescence (assuming a — 0.1, f E = 1, and e = 0). If the 
gap opens at larger separations, all the gas in the accretion 
disc around Mi is consumed before coalescence (assuming 
negligible refilling of gas through the gap). As reference, in 
the same figure we plot also a tr , to show the closest bi- 
nary separation at which we expect the gap to open. So, for 
Mi < 1 - 5 x 10 7 M e , assuming q = 0.01 — 0.1, the gap 
is opened well before t coa i = tduty, implying that the accre- 
tion is completed long before coalescence. Larger masses, in- 
stead, imply that gas accretion is still going on. In this case, 
a bright X-ray counterpart could be present in the error- 
box of LISA, given the large mass of the BHs. During the 
last year of inspiral, and during the BH plunge-in phase, 
gas dynamical perturbations and space-time curvature ef- 
fects strongly affect the disc structure. Armitage & Natara- 
jan (2002) have shown that, if some residual gas is present 
around the BHs just before coalescence, the accretion rate 
could be enhanced during the final GW driven inspiral. This 
process might trigger an episode of short-lived nuclear ac- 
tivity, so a transient, highly variable on-off precursor can 
anticipate the BH merger, making easier the identification 
of the galaxy hosting the LISA event. 

3.2 Tidal disruption of bound stars 

A burst of nuclear activity might be excited by the tidal 
disruption of a star bound to a BH, lasting ~ 1 yr (Rees 
1988). In order to produce an X-ray flare associated to a 
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Figure 1. The curve labeled with i coa j = tduty shows, as a func- 
tion of Mi, the largest BH binary separation for gap opening in 
order to have nuclear activity at coalescence. The curve labeled 
with t ooa i = tburst shows the separations below which the BH 
binary merges on a time scale tburst = 10 s yrs typical of a star- 
burst (see discussion in Sect. 4). The smallest separation at which 
the gap opens, ettn and a cr a are shown with dotted lines. Dashed 
lines correspond curves of constant coalescence time for GW emis- 
sion: 10 7 , 10 s , 10 9 , and 10 10 yr respectively, from bottom to top. 
We assume here q = 0.1, a = 0.1, = 1, and e = 0. 
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Figure 2. The same as Fig. 1 but for q = 0.01. 



LISA event, the star has to be disrupted not earlier than 
one-year before the binary coalesces. 

To asses the possible occurrence of such disruption 
within one-year prior coalescence two points must be sat- 
isfied: (a) the possibility for a star to survive in a orbit so 
bound that can interact with the secondary BH one year 
before coalescence; (b) the 3-body interaction must lead to 
the stellar tidal disruption instead of its ejection because of 
slingshot mechanism. 

To select binaries matching point (a), we can compare 
three length scales: (1) the tidal disruption radius of a so- 
lar type star orbiting around Mi rtd ~ 8.3 x 10" 6 pcMi'f; 
(2) the separation at which a massive BH binary coales- 
cence occurs in 1 yr: agw,bh ~ 3.7 x 10~ 5 pc Mi' 75 [q(l + 
q)] - 25 F(e) ' 25 ; (3) the separation at which a solar type 
star orbiting around the primary BH decays because of 
GW emission in, say, 10 7 yr which is the expected (short- 
est) lifetime of a BH binary in a galaxy nucleus: agw,* ~ 
3.7 x 10 pcMi®. The radius aaw,* gives the dimension 
of the empty-cusp around the primary BH. A star can be 
perturbed in its motion and be disrupted during the last- 
year of BH binary inspiral only if r t d < ^gw,bh ■ The re- 
gion in the Mi — q plane for which this condition holds is 
shown in fig. 3 as horizontally shaded area, and is bounded 
on the right by the condition rtd = 6GMi/c 2 : Note that 
for Mi > 2 x 10 s Mq a solar-type star has a tidal radius 
smaller than the radius of the last-stable circular orbit so 
that it would be swallowed by the BH rather than being 
disrupted. In addition, the existence of stars around the two 
BHs on the verge of merging, requires the size of the empty- 
cusp to be inside the BH binary separation at one-year from 
coalescence, i.e., qgw.bh > «gw,*. This region in the Mi —q 
plane is indicated with the vertically shaded area. The su- 
perposition of the two areas defines the parameter space of 
binaries that are allowed to host a star so bound to Mi that 
feels the presence of M2 during the last year of the binary 
inspiral. 

The point (b) is investigated by means of three body 
scattering experiments. A detailed description of the 3-body 
procedure can be found in Sesana, Haardt & Madau 2006. 
Here we briefly review those technical features which are 
of interest for the present application. The three masses 
are treated as point-like particles, and the tidal disrup- 
tion radius considered for the two BHs. We integrate the 
nine coupled, second order, differential equations of motion, 
assuming Newtonian gravity, and the mutual force acting 
on the two BH induced by gravitational wave in emission 
in quadrupole approximation (Iwasawa, Funato & Makino 
2005). We set the initial conditions as follows: we consider 
circular BH orbits, at initial separation <jgw,bh; bound stars 
are then drawn from an isotropic velocity distribution spher- 
ically symmetric around Mi, with the constrain that the 
distance from Mi at periastron is > r t d, and at apoastron 
is < cigw.bh — rtd,2, where r t d,2 is the tidal disruption ra- 
dius of M2. The integration is stopped as soon as one of 
the following conditions is met: (1) the star is kicked out by 
the binary with positive total energy, (2) the star crosses rtd 
or r t d,2- During the integration of the orbit, as the binary 
shrink because of GW emission, M2 perturbs the stellar or- 
bit around Mi. Typically the star is tidally disrupted by 
Mi or M2 well before it gets a positive energy because of 
slingshot mechanism. In a sample of about 5000 simulated 
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Figure 3. In the plane M\ — q, the region above the long dashed 
(short dashed, dotted-dashed) line limits LISA potential targets 
with SNR>5 at z = 1(3, 5). The horizontally dashed area marks 
systems with r tt j < agw,bh > the vertically dashed area marks 
systems with hgw.bh > ogw,* (see text for discussion). 



orbits, we recorded ~ 4950 disruption events, i.e. ~ 99% of 
the total. In other words, (almost) all binaries satisfying our 
condition (a) above, are also producing a tidal disruption 
event. 



3.3 Disc accretion: the afterglow in the era of 
next generation X— ray satellites 

Right after coalescence, the central relic BH is embed- 
ded in a hollow region surrounded by a gaseous disc. The 
gap will be re-filled on a time scale comparable to the 
gas viscous timescale. Milosavljevic & Phinney (2005) have 
shown that the AGN turns on after a time t on ~ 7(1 + 
z)(M B H/10 e Mq) 1 - 32 yr, creating an Eddington-limited X- 
ray source of luminosity L ~ 10 43 5 (M B h /10 6 Mq) erg s" 1 . 
Such a turn-on of nuclear activity could be the clear signa- 
ture of a very recent coalescence, and would allow for the 
unambiguous EM counterpart identification among other X- 
ray sources in the LISA error cube (see Hughes 2002, Vecchio 
2004 and Kocsis et al. 2005 for a discussion on the LISA spa- 
tial resolution). The next generation of X-ray missions (such 
a.s XEU& and Constellation- expected to be operating 
simultaneously to LISA, will be able to detect these objects 
up to 2 ~ 20 (Milosavljevic & Phinney 2005). 

The afterglow delay t on is a function of BH mass, and, 
as an example, only BHs lighter then 3 x 10 Mq can be 
observed during an operation time of ~ 30 yr. To quantify 



t www.rssd.esa.int/index.php?project=XEUS 
§ constellation.gsfc.nasa.gov 



the number of X-ray afterglows detectable during the op- 
eration lifetime of next generation satellites, we must rely 
on a model for the BH assembly in hierarchical models of 
structure formation. To this aim, we started with the results 
of Sesana et al. (2004, 2005), which calculate the rate of BH 
binary coalescences as a function of redshift and mass (see 
their fig. 1). The contribution of massive BH binaries to the 
LISA data stream was computed using an extended Press- 
Schechter merger tree code able to follow the assembly of 
massive BHs (Volonteri, Haardt & Madau 2003^) The to- 
tal number of coalescences that can be detected in a 3 years 
LISA mission assuming a signal to noise threshold for detec- 
tion S/N=5 is shown with a solid line in fig. 4, and turns out 
to be ~ 35, with a peak of detection around z ~ 5. Now, 
if we (i) assume that the disc spectrum is described by a 
thermal modified blackbody for a nearly maximal spinning 
BH, so that most of the luminosity is emitted at rest-frame 
energies hv ~ (0.5 — 5) keV (Milosavljevic & Phinney 2005); 
(ii) neglect photo-electric absorption in the X-ray observed 

band I , and (iii), take as fiducial value for future X-ray mis- 
sion a [0.5-2] keV flux limit ~ 10~ 1S erg s _1 cm~ 2 , we can 
compute the number of LISA sources detectable by XEUS 
in 1, 5, and 20 yrs operation time (fig. 4). Already after 1 yr 
from the LISA detection, almost 25% of the coalescence can 
be identified as X-ray sources. This fraction increases up to 
60% (80%) if longer lifetime for XEUS are considered. It is 
interesting to note that almost all the coalescences observed 
by LISA at z > 6 have an EM counterpart detectable by 
future X-ray missions within the first 5 yrs after the GW 
burst. At these redshifts, LISA error cube will be crowded 
by a large number of X-ray sources, so the sudden turn-on 
of a new X-ray source can give a peculiar fingerprint of the 
coalescence event. We also verify that increasing the flux 
limit (i.e., increasing the BH mass threshold for the after- 
glow identification) of future X-ray missions by a factor of, 
say, 10, our predicted counts remain valid, as afterglows are 
expected to have fluxes typically exceeding 10 -17 erg/s/cm 2 
in the X-ray band. 

As far as EM counterparts in other wavelengths are con- 
cerned, we can speculate that a wiggling radio jet could be 
observable during the last phases of spiral in of the BHs 
(see as a review Komossa 2006 and references therein) in 
particular when the precursor is present that guarantees the 
feeding of the BHs over most of their inspiral. This particu- 
lar feature could also persist for several years after the LISA 
detection. 

4 STARBURST GALAXY AS EM 
COUNTERPART 

In standard hierarchical scenario of structure formation, the 
assembly of massive BH is connected to galactic mergers. 

» Note that Volonteri et al. (2003) assume a short coalescence 
time invoking stellar dynamical processes with a continuos supply 
of low angular momentum stars. 

II This assumption is justified particularly for high redshift 
sources, since in this case we sample the high energy tail of the 
spectrum. 
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Figure 4. Number of coalescences as a function of redshift de- 
tectable by LISA in three years of operation (solid line; S/N 
> 5) and expected number of X-ray counterpart identifications by 
XEUS in 1, 5, and 20 yrs with dashed-dotted, short dashed, and 
long dashed line, respectively (neglecting absorption). The stan- 
dard ACDM cosmology with H = 70 km s _1 Mpc~\ Q m = 0.3, 
and f^A = 0.7 has been assumed. 
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Figure 5. Run of a/atr as function of evolution time. The vertical 
solid line marks a = a cr it/at r , while the horizontal line gives 
t burst . Different curves refer to different masses of the primary 
BH, labeled with log (Mi /M©). The top (bottom) panel shows the 
results for q = 0.1 (q = 0.01). We find that the starburst is likely 
to be still active at the coalescence only for Mi > 10 (2) X 10 6 Mq . 



Kazantzidis et al. (2005) have performed N-body SPH sim- 
ulations of binary equal- and unequal-mass (1:4) mergers of 
disc galaxies following the BH decay down to a separation 
< 100 pc, comparable to the numerical resolution. These 
simulations show an intense burst of star formation excited 
when the two BHs form a pair. Almost 90% of the central gas 
is converted into stars in less than tburst — 10 8 yr, implying a 
star formation rate of 30—100 Mq yr -1 , comparable to that 
observed in (ultra-)luminous infrared galaxies ((U-)LIRGs). 
The existence of such an intense burst of star formation as- 
sociated with a merger event (see also Di Matteo, Springel 
& Hernquist 2005; Springel, Di Matteo, & Hernquist 2005), 
could be, in principle, a peculiar feature of galaxies host- 
ing coalescing BH binaries. The peculiar correspondence is 
observable if BH coalescence is attained before starburst is 
completed, i.e., if BHs merge in less than 10 s yr. As seen 
in Section 2, dynamical friction against stellar and gaseous 
background is able to bring the BHs to the gap opening in 
only ~ 10 7 yr. As we already pointed out, for q < 0.1, the 
interaction between the BHs and the inner edge of the cir- 
cumbinary disc can lead to coalescence on a timescale i coa i 
given in Eq. (8). 



In Fig. 5 we show the run of a/atr with time t, assuming 
that dynamical friction have driven the BHs to a separation 
Oti (see eq. 3) in ~ 10 7 yr. The vertical solid line marks 
a = a cr it (eq. 5), the separation at which viscous torques be- 
come inefficient. For smaller separations, the evolution of the 
binary is dominated by GW emission. Different curves refer 
to different masses of the primary BH. The solid horizontal 
line indicates the typical starburst duration, tburst- The top 
(bottom) panel shows the results for q — 0.1 (q = 0.01). We 
find that the starburst is likely to be still active at the coa- 
lescence (i.e. £ coa i < tburst) only for Mi > 10(2) x 10 6 Mq. 
On the contrary, for smaller masses, the starburst is con- 
sumed before the BHs merge. In case dynamical friction is 
inefficient until a tI , the starburst activity at coalescence is 
still possible as shown in figs 1 and 2. The dot-dashed line 
marks, in the Mi — a plane, the locus t coa i = % urst i assuming 
that at these separations the gap is already open, and that 
the dynamics is dominated by viscous torques. Following 
this argument, the identification of a (U-)LIRG in the LISA 
error cube could be a probable EM counterpart of the GW 
event. Some caveat to this analysis has to be discussed: (i) 
the more massive BHs (> 10 8 Mq) are preferentially hosted 
in gas poor early type galaxies that could be the outcome 
of a gas-rich major merger. However there is the possibil- 
ity that dry mergers are important, i.e., mergers between 
gas poor galaxies, and in this case an intense burst of star 
formation is not expected; (ii) binaries with small q may cor- 
respond to small mass ratio of the interacting host galaxies, 
so that the associated starburst might be less powerful then 
what typically observed in (U-)LIRG; (iii) the presence of 
an intense starburst event in the last phases of galaxy merg- 
ers could depend on the morphology of the two galaxies and 
in particular on the relative size of the bulges, poorly ad- 
dressed by numerical simulations up to now. Nevertheless, 
the presence of a starburst galaxy in the LISA error cube 
may indicate a preferential host of the BH coalescence. 
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Figure 6. Summary sketch of the main results. Nuclear activ- 
ity as a precursor is expected to occur in binaries with primaries 
heavier than ~ 10 7 Mq . Coalescences involving less massive BHs 
will lack of any precursor and, for Mi < 3 X 10 6 Mq, are followed 
within 20 yr by an EM afterglow. The small triangle labeled as 
"tidal" shows the region of the parameter space where tidal dis- 
ruption events of bound stars may be associated to the LISA 
detection. Finally, an intense starburst event may be still present 
at the coalescence of massive BHs (Mi >2x 10 6 Mq) embedded 
in gas rich galaxies. 



Smaller BHs (Mi < 3 x 10 6 Mq) will lack of any pre- 
cursor, since we have found that the gas around the BHs is 
consumed long before coalescence. Moreover, the starburst 
activity is also unimportant, since the merging time of the 
BH binary is longer than the typical starburst timescale. 
By contrast, a LISA event involving light BHs is followed 
by an afterglow with rise-time less than 20 yr, i.e. the ex- 
pected operation time of XEUS. Using cosmological merger 
tree algorithm, we have shown that XEUS will be able to 
identify almost all the predicted LISA events in 20 years of 
operation, and the totality of the events at z > 6 in only 
5 years. There is a small transition region in between the 
light and heavy BHs considered in this paper, where a LISA 
event may be only associated to a starburst but shows no 
precursor. 

The association of an EM counterpart to a LISA event 
should be linked to a variable on-ofT (precursor) or off-on 
(afterglow) state. In addition, we predict that the precur- 
sor is related to mergers of heavy BHs, whereas afterglows 
are detectable only for lighter BHs. This dichotomy may be 
tested by comparing the redshift of observed EM counter- 
part hosts with the redshifted masses estimated analysing 
the LISA data stream (see, e.g., Vecchio 2004). 
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5 SUMMARY AND CONCLUSIONS 

We have explored a number of possible EM counterparts of 
the coalescence of massive BH binaries at cosmological dis- 
tances, assessing their detectability with future X-ray mis- 
sions. Our main findings are summarized in fig. 6. 

EM precursor is expected to occur in binaries with pri- 
maries heavier than ~ 10 7 Mq. If this event is related to a 
merger of quasi equal mass, gas rich galaxies, star formation 
activity should be still present at the coalescence time. In 
this case, the identification of a starburst galaxy in the LISA 
error cube could be a probable EM counterpart of the GW 
event. 

AGN precursors, detectable during the last year of in- 
spiral, can be identified as variable X-ray sources switch- 
ing to an off state. X-ray flares due to accretion of matter 
pushed by the incoming secondary BH might be the finger- 
print of the LISA event. This on-off preglow would help in 
identifying the host galaxy within the error cube. Note that 
the extrapolation of our results to q > 0.1 are only specula- 
tive, and demands more detailed studies (solid lines ending 
with dots in fig. 6 indicate uncertain extrapolation) . EM af- 
terglows associated to these massive binaries require onset 
times in excess of > 100 yr, and so exceedingly long expecta- 
tion times. The triangle involving the heaviest BHs denotes 
the region where a tidally disrupted star may be accreted in 
the last year before coalescence, potentially causing a X-ray 
episode of accretion. This would again appear as a preglow. 
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